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Abstract 

This paper presents and compares three analytical methods for calculating low 
frequency band gap boundaries in doubly periodic arrays of resonating thin elastic 
shells. It is shown that both lattice sum expansions in the vicinity of its poles and 
self consistent scheme could be used to predict boundaries of low-frequency band gaps 
due to axysimmetric (n = 0) and dipolar (n = 1) shell resonances. The accuracy of 
the former method is limited to low filling fraction arrays, however the application 
of the matched asymptotic expansions could significantly improve approximations of 
the upper band gap boundary due to axysimmetric resonance. The self-consistent 
scheme is shown to be very robust and gives reliable results even for dense arrays with 
filling fraction higher than 40%. The results are used to predict the dependence of the 
position and width of the low frequency band gap on the properties of shells and their 
periodic arrays. 

1 Introduction 

Periodic and random arrangements of scatterers with low frequency resonances in air have 
attracted a significant amount of attention recently due to their applications as locally 
resonant acoustic metamaterials [H [21 [3]. Different types of scatterers have been shown 
to possess these types of resonances, among them most notable are split rings i.e. 2D 
analogues of Helmholtz resonators [I] . 

It has been demonstrated recently that periodic arrays of thin elastic shells in air can 
possess multiple resonant band gaps in the frequency range below the first Bragg band gap 
[U [6]. Data from experiments with finite arrays of elastic shells [5] have shown that the 
insertion loss (IL) peak corresponding to the axisymmetric (n = 0) resonance is comparable 
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to that for the first Bragg band gap but at a lower frequency. The IL peak corresponding 
to the n = 1 resonance also appears at a lower frequency but is less strong. 

The position and width of the complete band gaps can be obtained using standard nu- 
merical and semi-analytical techniques [3(8]. However it would be attractive to have analyt- 
ical or semi-analytical expressions allowing approximate estimates of band gap boundaries. 
This would help to better understand their dependence on scatterer and array parameters 
and hence facilitate the design of resonant metamaterials with desired properties. 

In this paper several methods allowing simple estimations of resonant band gaps bound- 
aries in periodic arrays of elastic shells are presented. In Section [2] the dispersion relation 
and expressions for lower and upper band gap bounds are derived using expansion of the 
lattice sum within the vicinity of its pole [9]. The derived dispersion relation is similar to 
that obtained with the Foldy's approximations [12]. This method allows analytical approx- 
imations of two resonant band gaps (n = and n = 1). In Section [3] matched asymptotic 
expansions are used to improve approximations for the upper boundary of a band gap 
associated with the axysimmetric (n = 0) resonance. A self-consistent method, also known 
as coherent potential approximation (CPA) |10(. Ill j . is used to derive simple expressions for 
effective density and bulk modulus of the array and hence to obtain a dispersion relation 
in Section HJ The band gaps due to n = and n = 1 resonances are determined by finding 
frequency ranges where effective density or bulk modulus (respectively) are negative. The 
results of all techniques are compared with each other and with numerical model predic- 
tions. In Section [5] the derived approximations are applied to estimate the width of the 
band gaps with respect to the variation in elastic shell radius, its thickness, lattice cell size 
(i.e. lattice constant) and Young's modulus of the elastic shell. 

2 Foldy's type formula 

First consider an eigenvalue problem stated for acoustic wave propagation through a doubly 
periodic array of circular scatterers Cj, j = 0,1,2.... The scatterers are arranged in an 
infinite lattice A with lattice constant L. The acoustic environment is characterised by 
its density p Q and sound speed c Q . In this section the origins of Cartesian coordinates 
r = (x, y) and polar coordinates coincide with the center of scatterer Cq in the primary 
cell. Centre Oj of each scatterer Cj in the jth cell of lattice A is defined by the position 
vector 



where a\ and a? are the fundamental translation vectors. The solution in acoustic medium 
p(r) satisfies the Helmholtz equation, 



Rj = n\a\ + n2d2, rai, ri2 G Z, 
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where A = — — I r— H — ^tt^t, k a = lo/c is the wavenumber defined as a ratio be- 
r or \ or J r o9 2 

tween angular frequency uj and sound speed c . Throughout this paper the time-harmonic 
dependence is taken as exp(-iujt). 

Solution p(r) is subject to the boundary conditions specified later in this paper and 
the Floquet-Bloch conditions, also referred to as quasi-periodicity conditions, that are 

p(r + Rj) = p(r); (3) 

where r is the position vector of a field point, Rj is the position vector of each scatterer 
centre and (3 = {qi,q 2 ) is a given wave vector. The components of the wave vector f3 can 
also be defined in polar coordinates as q\ = f3 cost and q 2 = (3smr. 
The general solution of the Helmholtz equation ([2]) is given by 

+oo 

p(r,6)= Y, [B J n Jn{k r) + BlY n {k r)}e ine . (4) 

n=—oo 

Applying the quasi-periodicity conditions (J3]) together with the boundary conditions im- 
posed on the surface of each scatterer to equation Q results in the dispersion relation 
referred to as the Rayleigh Identity j7j 

B n ~ E (-l) p - n Vp-n(ko,P)Z p BJ = 0, n G Z (5) 

p= — OD 

where factor Z n is defined by boundary conditions on scatterer's surface. For example, 
zero normal velocity at the surface of the rigid scatterers of radius a gives 

J' n (k a) 

Zn -~Y>JM () 

with prime denoting derivative with respect to polar coordinate r. The lattice sum 
a^(k ,f3) in equation ([5]) can be derived in the lattice A as 

a^(k D ,(3)= £ e^Y n (kR )e ina ^ (7) 
RjeA\{0} 

where representation Rj = Rj (cos otj , sin aj ) has been used. 

The lattice sum can also be derived in a reciprocal lattice A* defined by lattice vectors 

R* m = 27r(mi&i + m 2 b 2 ), m = (mi, m 2 ) G Z 2 . (8) 

The fundamental translational vectors b± and &2 m A* are related to those in A through 

a,ibj = dij, i,j = 1,2. (9) 
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It is shown [13] that 



where A is the area of a lattice cell,/3 m = (3 + = /3 m (cosr m ,sinr m ) and £ G [0, £] with 

C ^ min R j. (11) 

The lattice sum given by (fTUj) has simple poles /c = /3 m that correspond to a plane 
wave solutions satisfying Floquet-Bloch conditions ([3]) in the absence of the scatterers [9]. 
The magnitude f3 m is defined by M pairs of integer numbers (777,1,777,2) G T? M that gives M 
different wave vectors (3 m . In the vicinity of this pole the related M terms of the lattice 
sum (fT0|) take the leading order that is 

o*M) = ^ Y, e ~W^ 1 + ^ , ( 3), (12) 

A K o Pm 



m£Z 2 



M 



where On ,2 (k 1 (3) is the next order in expansion over 

(k 2 - p^n)L 2 < 1, me 1 2 M . (13) 

In this section it is assumed that only one pair (7711,7712) = (0,0) contributes to the 
leading order of expansion (|12p . This results in (3 m = f3. It is noted that this is true for 
the case when k Q L <C 1 and f3 m L <C 1. 

It is also assumed that j3L is taken within the first irreducible Brillouin zone [8j that 
can be described by the contour (0, 0) — (tt, 0) — (ir, tt). 

The smallness of (3L and approximation of lattice sum (|12p by the summand related 
to the pair of (mi, 7712) = (0, 0) make valid further approximations along (0, 0) — (ir, 0) and 
(tt, tt) — (0, 0) lines. This also leads to a symmetry around point (0, 0) so that vector f3 is 
only considered within the interval (0, 0) — (tt, 0) (i.e. (3 = q\ and r = with q\ G [0, q] 
and qL <^ ir). As a result the leading order of the lattice sum ([12]) is simplified to 

i n 4 

a n(.k,f3) « a^' 1 (k,(3) = -j-j^r^, ( 14 ) 

Substituting equations (|14p into dispersion relation ([5]) one can obtain the following 
system of homogeneous algebraic equations with respect to unknown coefficients B^, n G Z 

4 +°° 

B l ~ A(k 2 _ p) £ ^Y~ n B J p Z P = 0, n G Z. (15) 
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By multiplying both sides of equation (fT5|) by (— i) n it is possible to redefine the unknown 
coefficients as = (—i) n B^, n G Z so that the system of equations (|15p becomes similar 
to that derived in [151 ec L- (71)]. This equation can now be employed to derive Foldy's type 
dispersion relation. Following the derivations in \15\ eqs. (76) and (77)] it is concluded 
that coefficients B^, n G Z are identical for all indices n G Z, yielding 

/? 2 = fc' " (16) 

+oo 

where F = Z p is determined by boundary conditions on the scatterer's surface. It 

p=— oo 

also coincides with the imaginary part of the far-field pattern of a single scatterer. The 
consistency of equations (|13p and (fTUj) requires that 

^-FL 2 « 1. (17) 

It must be noted that for the resonant scatterers F have singular points that give reso- 
nance frequencies. Therefore, in the vicinity of resonances the asymptotic orders of factors 
Z n , n G Z can invalidate the derivation of dispersion relation (]16p . This, however, can also 
signify the existence of a band gap due to resonance. 

Now consider a square periodic array of thin-walled elastic shells with area of lattice 
cell A = L 2 . The application of long- wave low- frequency approximations results in the 
following expressions for factor Z n , n G Z described by [U eq. (22)] 

7 = \J'n(koa)?[l - (ha) 2 + n 2 ] 

J4(A: a)^(A :o a)[l-(fc3a) 2 + n2] + [n2-(A ; 3a) 2 ] /9o ( /9 7ra/ i )- 1 ' ~ n n ' 1 ; 

where p is the density of an elastic material, £3 = w/c3 and a is shell mid-surface radius 
(referred to as radius of the elastic shell) and h is its half-thickness. Here the dilatational 
wave speed C3 for a thin elastic plate is recalled 



E 



within which S and z/ are Young's modulus and Poisson ratio of elastic shell, respectively. 

For the low frequency range that contains n = (i.e. axisymmetric resonance) and 
n = ±1 resonances the far-field pattern of elastic shell can be approximated by its first two 
terms that gives 

F = Z + 2Z 1 . (20) 

The assumption that wavelength of a propagating acoustic wave is much bigger than the 
lattice cell size (i.e. lattice constant L) immediately results in 

e = k a < 1. (21) 
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The latter can be used to simplify (|T8|) by expanding Bessel functions. It is also assumed 
that the elastic shell is "soft" (i.e. p Q /p = 0(e 2 ) and c^/c = 0(e)) so that acoustic waves 
are able to penetrate it. Thus 



Zn 



2*K 2 -k 2 



»7T 



k 2 



p Q a Kq - k 2 



o +£ 2 



where wavenumbers 



ph 2K 2 - k 2 



a 

Kq = -y/K-p ~ 



^ log ^ + ^(5 + 4 7 ) 



(22) 
(23) 

(24) 
(25) 



with 



and 



Po_a 
p h 



(26) 



(27) 

correspond to axisymmetric resonances of shell in vacuum with vacuum and acoustic 
medium inside, respectively. In equation (j23|) approximated Z\ has a singularity at k 2 = 
2Kq related to the leading order of resonance frequency with index n = 1. From this 
approximation it can be concluded that frequency of n = 1 resonance of a thin elastic shell 
defined from 

K x = V2K (28) 

has no dependence on the physical parameters of the surrounding acoustic environment in 
the leading order. 

As f3L — > equation (|16[) gives two non-zero solutions for k Q in addition to k a = 0. 
These solutions can be used to find the size of the corresponding band gaps. 

In the vicinity of axisymmetric resonance Kq as {3L — > the dispersion relation (|16p 
can be transformed to 

TK P (K C - K p ) + [(k a) 2 - (JC C + K p j\ [(1 + F){K P - JC C ) - TJC P {1 + 2JC P )] = 0. (29) 

This results in the following solution for k D > 



«o,2 



'/C p + /C c + — 



TJCpQCp - JC C 



a 2 K 2 + FK P 



(30) 
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Figure 1: (Color online) Foldy's approximation (116[> ('- -') compared with the semi- 
analytical solution of Rayleigh Identity ([5]) (solid line), (a) a = 0.0275 m and (c) a = 0.0375 
m. 



where T = no 2 /A is the scatterer filling fraction. Expression (|30p immediately estimates 
the upper limit of the band gap. Note that the lower limit of the mentioned band gap 
is given by fc 0i i = KqL. Therefore, based on the derived approximations, it is possible 
to conclude that in the doubly periodic array of thin elastic shells (with filling fraction 
T <C 1) waves do not propagate in the frequency interval 



Jn=0i J n=0 

In the vicinity of n = 1 resonance K\ is approximated by (|28p and as PL — > the 
dispersion relation (]16p can be transformed to 

2FK p K, c {K, p - K c ) + [{k af - 2/C c )] [(1 + 2F- 2FK P ){1C C - K p ) - FK C {1 + 2K P )\ = 0. 

(32) 

For k Q > the solution of this equation takes the following form 




2J-K, p K, c (K, c — KL p ) 
K c -K p + T [2(/C c - K P ){1 - Kp) - JC c (l + 2)C P )} 



« ^a 2 Kl + 2FK p lC c (33) 

The lower limit of the band gap related to n = 1 resonance is given by equation (|28|) . 
Therefore the frequency interval where waves do not propagate is estimated as 



Jn=l ) In- 



_C3_ 

2na 



1 + 



1C P c 3 



1 + ^(1+ JO 



(31) 



\/2c 3 y/2, 



C3 



27ra 27ra 



(34) 
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In FigureQ]the solutions of dispersion relation (|16p obtained for F that contains only two 
terms n = 0, 1 are shown. Throughout this paper unless otherwise specified the material 
parameters of thin elastic shell are identical to those in [5] and its thickness 2h = 0.00025 
m. The radius of the elastic shell is varied whereas the lattice constant L is fixed (L = 0.08 
m) so that two values of the filling fraction J- are considered: 0.4 and 0.6. Two band gaps 
related to the resonances Ko, n = and y/2Ko, n = 1 are observed below the first band gap 
associated with the array periodicity. According to Figure QJa) the estimates of the lower 
and upper limits of the band gap are within 5% of the exact values as long as the filling 
fraction is relatively low (J 7 > 0.4). However as the filling fraction is increased up to 0.6 
the accuracy is expected to decrease according to the condition (fTTl) . This is illustrated in 
Figure DJb) where estimated value of the upper boundary of the band gap associated with 
n = is 30% in error from the exact value. It must also be noted that this upper boundary 
obtained with Foldy's equation depends on the form of the factor Z n . Its approximated 
form might give more accurate results than the original form of Z n given by (|18p . This 
inconsistency points out that to obtain more accurate and consistent results equation (1161) 
has to be modified to include the periodicity effects. 

3 Matched asymptotic expansion 

The results obtained in the previous section can be improved by using the technique based 
on matched asymptotic expansions (MAE) [91 [16] . This technique allows approximating 
eigenvalues of ©-([3]) when k Q L = 0(1) in the vicinity of band gaps associated with the 
periodicity (Bragg band gaps). However in this paper MAE is used to estimate the solution 
that forms the upper boundary of the band gap due to the axisymmetric resonance (n = 0) 
of a thin elastic shell. 

First outer and inner regions are introduced. In the inner region surrounding each 
scatterer the solution is solved in conjunction with the boundary conditions imposed on 
the surface of the scatterer. The characteristic length of this region is the radius a of 
scatterer. Thus a dimensionless inner coordinate is given by 



In the outer region the scatterers are replaced by point sources and the solution is subject 
to quasi-periodicity conditions. The characteristic length in this region is equal to the 
wavelength of sound. Thus a dimensionless outer coordinate is given by 



The wavelength is assumed to be much bigger than the scatterer radius. This requirement 
results in small parameter e introduced in Section [2] that connects two regions through 



r 

x = — . 

a 



(35) 



y = k Q r. 



(36) 



y = ex. 



(37) 
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3.1 Inner solution 

The inner solution is found from the boundary value problem for a single thin elastic shell. 
The solution around the cylinder is represented as 

+00 

il>(x,9)= B J n [J n {ex) + Z n Y n {ex)\e me (38) 

n=— 00 

with unknown constants and coefficients Z n defined by equation (|18p . According to 
equation (f22|) . factor Zq takes order 0(e 2 ) if its singular points (i.e. resonances) are isolated. 
In the vicinity of axisymmetric resonance the order of factor Zq can be transformed to 0{rf) 
that is 

Z = 5 ZV , with 8 z = 0(l), (39) 
where small parameter 77 is bigger than e and is assumed to be 

V = — ^ , (40) 

K — log e 

within which K is the unknown constant. 

By assuming Zq of order r/ the proximity of ko to the axisymmetric resonance Kq is 
restricted by 

K - Ko = O (£) (41) 

Expanding solution (|38p up to the order 77 first with respect to the inner coordinate x 
and then with respect to the outer coordinate y results in 



^)(x,9) = B J Q 



2 2 
1 5z + ri—5z logx 

7T 7T 



(42) 



It is noted that the terms in O(l) and 0{rf) orders have no azimuthal dependence. 



3.2 Outer solution 

The outer solution of the Helmholtz equation ([2]) satisfies the quasi-periodicity conditions 
© and is singular at the points Oj. This results in 

+00 

tf(r,0)= A " E e^H^ikorj)^, (43) 

n=-oo Hj-eA\{0} 

where local coordinates (rj,6j) with origin at Oj are used. 

It is assumed that the leading order of coefficients A n , n G Z is bounded by the small 
parameter 77. This approach is similar to that derived in p3] and |16t Section 6.8]. Thus 

A n = rjA n . (44) 
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Figure 2: (Color online) MAE approximation of the upper limit (I55jl ('• • •') compared with 
the Foldy's approximation (|16|) (' — ') and the semi-analytical solution of Rayleigh Identity 
© (solid line), (a) a = 0.0275 and (b) a = 0.0375. 



The use of addition theorem for the Bessel functions [17] and lattice sum representation 
in the reciprocal lattice A* transform the outer solution (|43p to 



+oo 



+oo 



*W(y,0)=ri £ A n \HW(y)e me + £ (-l) n -^ p (k , (3)J p (y)e^ e 



(45) 



p= — OD 



and 



cr n (k ,(3) = -5 ,n + ia%(k ,J3), (46) 

where a^(k ,/3) is given by equation (fTUl) . 

The form of the inner solution (|42p assumes that the leading order should include 
axisymmetric terms only. In contrast to Section [2] it is assumed here that for n = the 
entire lattice sum (fTO]) (not only pair (mi, 7712) = (0,0)) contributes to the outer solution 
as 

a%(k ,P) = -, with 6 = 0(1). (47) 
77 

Then expanding outer solution in terms of the inner coordinate x up to order 77 leads to 



2 2 

5 h 77— (K + 7 — log 2 + log x) 

7T 7T 



(48) 
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3.3 Matching 



The inner expansion of the outer solution (|48p can now be matched with the outer expansion 
of the inner solution (|42h . First factors of log £ in order rj give 

(49) 
(50) 
(51) 



A = -iS z B^. 
Matching constants in order r\ results in 

K = log 2 - 7. 

Finally the consistency of the leading orders in (|48p and (|42p requires that 

Z a%(k o ,P)-l = 0, 

where equations (f3"9"j) . (|17|) and (|4"9"j) have been used. 

As /3L —7- the lattice sum in equation (|51[) can be approximated by 



<tf(A: O) /3L->0) 



1 



Jo(U) 



4L 2 1 4L 2 x ^ 

A (A; L)2 o( ^ j+ A (k L) 2 -(R* m L) 2 



(52) 

It is noted that the last term given as a series over the reciprocal lattice represents the 
correction factor to the upper limit of the band gap related to the axisymmetric resonance. 
This correction factor includes the periodicity effects. 

For the square lattice with lattice constant L and £ = L/2 the lattice sum ([52]) is 
transformed into 



<r$(ko,PL = 0) 



1 



HKL/2) 



{KLf 



Y (k o L/2) + 4S(k a L) 



(53) 



with 



S(k L) 



J (vrymf+mf) 



(k a L) 2 -Air 2 {ml + 



ITln 



(54) 



m^(0,0) 



It is noted that £ is chosen from the interval set in (|lip to improve the convergence of S. 

The lattice sum ([53]) and the dispersion relation (f5"Tj) can be used to improve the upper 
limit of the band gap found in (f3"Uj) . Thus the improved estimate is given by the solution 
of the following equation 



Z {4 - (k L) 2 [Y (k o L) + AS{k L)]} - (k o L) 2 J (k o L) = 0, 



(55) 
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within which k Q L ^ and Jo(k L) ^ 0. To find the solution of equation (|55|) the infinite 
sum (|54|) has to be truncated. The results are accurate to three decimal places for both 
\mi\ > 8 and \m2\ > 8. 

Figure [2] demonstrates the improvement in approximation of the upper boundary of 
the band gap due to axisymmetric resonance given by equation ([55]) compared to Foldy's 
approximation equation (130 j) . For a low filling fraction (i.e. J- < 0.4) the results are 
accurate within 1% from those obtained with the Foldy's equation, see Figure (2)(a). As 
the filling fraction increases the approximation of the branch representing the upper bound 
deteriorates. However the limiting point at f3L = and its vicinity are not affected, see 
Figure [2Kb) so that the estimate remains valid even for large filling fraction (J 7 > 0.4) with 
the accuracy within 1% from the exact solution. 

4 Self-consistent effective medium formulation 

The effective medium approach provides an alternative method for estimating the bound- 
aries of band gaps due to resonances with indices n = and n = 1 . 

Consider an effective medium as a replacement of a periodic array of scatterers in a 
fluid matrix (see Figure [3]). The homogenization |10] 118] becomes possible only when the 
wavelength in both fluid and effective medium is much bigger than the radius of scatterers a 
and separation distance between them (lattice constant L) that is k a L <C 1 and k e ffL <C 1. 

The solution of the Helmholtz equation 

Ap a + k 2 aPa = 0, (56) 

is given in terms of the displacement potential p a (r) where k a = co/c a is the wavenumber, 
c a = v B a / p a is speed of sound, index a relates solution p to one of the regions (i.e. 'e//' 
effective medium (I), 'o' matching fluid layer between scatterer and effective medium). 
The problem is solved in polar coordinates r = (r,9). It is also noted that the physical 
parameters of the matching fluid layer coincide with the acoustic medium introduced in 
Section [2] 

To find the effective medium parameters (effective density p e ff and bulk modulus -B e //) 
and hence its wavenumber k e jf = ujy/p e ff/B e ff the problem is split into three regions that 
are [E] 

Region (I) - circular scatterer, 

Region (II) - fluid layer with properties identical to those of the fluid matrix, 

Region (III) - effective medium with yet unknown properties, 

where regions (I) and (II) are introduced to derive parameters of the effective medium 
through matching acoustic potentials at the interface of the region (II) and (III) r = R a . 
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Figure 3: Geometry of the effective medium 

That is 

Po = ^Peff> 

dPo = d^j_ (57) 
or or 

where outer radius R Q is derived in terms of the scatterer radius and the filling fraction of 
the original periodic array with lattice constant L so that the filling fraction in a composite 
inclusion remains equal to that of the original periodic array. This gives 

a 



R 



The total wave field in region (III) is given by 

Peff = Pine + P scat ) 

„. _ ik eff rcos(9~l3) 
Pine — c j 



(58) 

(59) 
(60) 



where pi nc is plane wave incident at an angle (3 and p sca t represents scattered wave field. 
Having assumed that the effective medium behaves as a homogeneous fluid the scattered 
wave field vanishes that leads to 

Peff = Pine (61) 

with the plane wave expanded over the regular Bessel functions J n (k e ffr) as [T71 eqs. 9.1.44 
and 9.1.45] 



Pine = ^2 inj n(k e ffr)e 



in(6-/3) 



(62) 
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Solution inside the annular layer (II) takes the following form 



+00 



Po= Yl B n [J n (k r) + Z n Y n (k r)]e ine , 



(63) 



where factors Z n are defined by (fT8|) , 

Substitution of equations ([62]) and ([63]) into the boundary conditions ([57|) leads to 



B n [Jn(koRo) + Z n Y n (k R )} = ^Hj n (k eff R )e in ^ 2 -^, 

Po 

Bn \JnQ^oRo) ~\~ Z n Y n (k R )] = J n {k e f fR a )e ^ ^\ 



(64) 



where n 6 Z. To solve system (|64p it first has to be truncated at some integer number iV 
that gives 



Ax 



(65) 



where matrix A has 4iV rows and 2N columns, vector x has 2N elements and vector b has 
4iV elements, and 



/ [J- N (koRo) + Z-NY- N (koRo)] 

[j'_ N (k R ) + Z- N YL N (k Ro)] 



\ 



V 

/ B^ N \ 

B-N+l 








[Jjv(fc i2o) + ^^V(^o^o)] 

[J^(fe i? ) + Z N Y^(k R )} J 



(66) 



a; 



(67) 



V Bn J 

( *££j_ 



ii )e-^^ 2 -« \ 



J^(A; e//J R )e-^^ 2 -« 

Mlj N {k e}f R )e iN ^^ 
V J^/^)^ 2 ^ 



(68) 



/ 



According to Kronecker-Capelli theorem [2T] the overdetermined system (|65p is com- 
patible if the rank of its coefficient matrix is equal to that of the augmented matrix. On 
the other hand the unique solution of system (|65p exists if the rank of matrix A is equal 
to number of variables (i.e. number of columns 2N). Hence combination of these two 
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statements and application of Gaussian elimination algorithm to the augmented matrix 
(A|6) yield the criteria for the existence of the solution that is 



J' n (k R ) + Z n Y^(k R ) _ p J' n {k e ffR ) 



n 



J n (k R ) + Z n Y n (k R ) p ef f J n (k e ffR ) ' 
These equations can now be simplified by recalling the assumptions 

7] = k R <C 1 and k e ffR a <C 1. 
Also equation (j69j) needs to be rewritten as 



■N..N. 



J n (rj) + Z n Y n (r]) j n {v£c) 



J n {r]) + Z n Y n (r]) J n (j)£c) 



n 



-N..N, 



(69) 



(70) 



(71) 



where '•' stands for the derivative of Bessel functions with respect to dimensionless param- 
eter rj. Parameters £ p and £ c depend on the effective medium properties as 



Po 

Peff'' 
Co 

°eff' 



(72) 
(73) 



Bessel functions and their derivatives in ()7ip can be replaced by their approximations in 
the leading orders as rj — > J7] , yielding 



ri\ n 1 

o) "7' n - ' 
2/ nl 
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— n 
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2 




nl 
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v -n-1 


2 







7T 



2tt 



-, n > 0, 
1 

2(n + l)! 7 ' 
-, n > 0. 



n > 0, 



(74) 
(75) 
(76) 

(77) 



It is convenient to consider subsequently that Bessel functions have only zero and positive 
integer orders with lo( 7 ?) = 2/7rlog(ry). For the negative orders of Bessel functions f n (rj) = 
{ J n (rj), Y n (r])} equation ([7T]) is identical to that with n > due to the relation f- n (i]) = 

Substitution of equations (j73]l - (177]l into equations (J7T]l gives for n = 



(ry/2) 2 - Zq/tt = t t2 (V^ 2 
l + 2Z Iog(»7)/7r ^ C V2, 



(78) 
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and for n > 



( V /2) 2n [(n - l)!]- 1 {1/2 - (rj/2) 2 [n(n + l)]' 1 } + Z n n\(2^)~ l _ g p 
(?7/2) 2 »*(n!)- 1 - Z n (n - l)!^ 1 2 



n 



2 7 n+l 

(79) 



Following the orders of smallness involved in equations (|78|) and (|79p the factor Z n has to 
appear as 

f (»7/2) 2 irZo,Q for n = 0, 

Z «= ( V /2) 2n \ Z nfi for n>0. ( 80 ) 
L (n — ly.nl 

This assumption can be easily justified for the circular scatterers by expanding Z n with 
respect to the small parameter, see ([22]) and ([23|h 

By collecting the same orders of smallness in (|78j) and (|79|) the leading orders can be 
derived in the following form 

^1 = 1 - Z„,o, n = 0, (81) 

^ = Jj^f^ - > 0. (82) 

From these two equations the parameters of the effective medium are derived in the fol- 
lowing order 

B eff = B 1 , n = 0, (83) 
1 — ^0,0 

Pe// = Po ! , f"' . n>0, (84) 

1 + ^n,0 

where B Q = c 2 p Q is bulk modulus of matching fluid layer (i.e. air). 

The effective density (I84D depends on index n since coefficient Z n $ ~ T n . This limits 
the use of the proposed model and can only include the contribution of two harmonics 
that are n = and any non-zero index n. The usual choice of n = and n = 1 which 
are the harmonics that contribute most is followed now. Using (|22p and (|23p the following 
expressions are derived for effective density and bulk modulus: 

B e ff _ 1C C + 1C P - (k Q a) 2 ^ 



B JC c {l-T) + lC p -{k a) 2 {l-Ty 
p e ff _ K c [2 + F{2 - K p )} - {k a) 2 [l + F(l - K p )] 
Po K c [2 - F{2 - K. p )} - (k a) 2 [l - T(l - K p )] ' 

where JC P and JC C are defined by (f26j) and ((3j 



(86) 
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Figure 4: (Color online) Dispersion relation predicted by self-consistent method (187p ('- • - 
') compared with the semi-analytical solution of Rayleigh Identity ([5]) (solid line), Foldy's 
approximations (US) ('- -') and MAE solution ([55]) ('• • •')• (a) a = 0.0275 and (b) a = 
0.0375. 



The dispersion relation of the effective medium is 



k 



eff 



Peff B o r,2 



p (87) 

which has two poles corresponding to resonances with n = and n = 1. Moreover there 
are two intervals where k e ff is imaginary, i.e. band gaps. In the first interval the bulk 
modulus B e ff in (|85j) is negative, 



1 



k Q > -^/K,p + K, c = K , 
a 



a V 1 - T a 



pi 



whereas the second interval gives negative effective density p e ff, 



k a > - J 2K, C 



1 



K < 



'2/C c + 



1 - F(l - K,p) jr <x a 
1 



-\la 2 K\-TKpK c 



FK,pK, c 



1 + F{1 - K p ) t«i a 



i 2 K\ + F)CplC c . 



(88) 
(89) 

(90) 
(91) 



Therefore the band gap due to the axisymmetric resonance corresponds to the frequency 
range where effective bulk modulus is negative: 



Jn=0i Jn=0 



_C3_ 

2-Ka 



1 + 



K, r 2ira 



1+ l 



1 + 



F 



l-F 



(92) 
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while band gap due to n=l resonance corresponds to the frequency range where effective 
density is negative: 



J 71=1 ) J 71=1 



V2c 3 J T1C P V2c 3 

2na V 2 [1 - T(l - K p )\ ' 2vra ' 



1 + 



2[1 + ^(1-/C P )] 



(93) 



Comparing n = band gap limits approximated by (|92p with those derived using lattice 
sum expansion (I3ip it can be seen that they coincide as J- — > 0. On the other hand the 
limits (|93p for the n = 1 band gap are shifted towards lower frequency compared to the 
interval (134p . It is noted that as J- — > the approximation of the lower limit f n= \ coincides 
with the resonance n = 1 of the elastic shell embedded into the hollow cylinder with rigid 
walls and radius equal to that of the composite inclusion. 

In Figure 4 the dispersion relation (|87p is compared with the approximations derived in 
Sections [2] and [3] in the interval k e tfL £ [0, tt]. This interval coincides with that chosen in 
Section [2] for the wave vector (3 (i.e. (3L belongs to [0, n]). It is observed that for relatively 
low filling fractions (see Figure IH(a) for T ~ 0.4) as well as for high filling fractions (see 
Figure E^b) for T ~ 0.6) the results are accurate to within 5% of the semi-analytical 
solution. It is also noted that as the filling fraction increases the approximation of the 
upper limit of n = band gap is less accurate than that obtained with equation (I55p . 



5 Results. Band gap width 

In Figure [5] the results illustrate the application of the approximations derived in the 
previous sections to estimate the width of the band gaps. 

The dependence of the band gaps related to the shell resonances with indices n = 
and n = 1 on the filling fraction is shown in Figures [5]^ a) and (b). The Bragg frequency 
/ = c /(2L) is constant with respect to T in Figure [5{a) where the characteristic size of 
the lattice cell (i.e lattice constant L) is fixed and the shell radius is varied. It is noted that 
the approximations based on the Foldy's equation underestimate the upper limit of n = 
band gap as J- — > 1. This is improved with the help of the matched asymptotic results and 
the self-consistent method. It is also observed that for low filling fractions (i.e. J- < 0.1) 
the axisymmetric resonance and the corresponding band gap are expected to be observed 
above the first Bragg frequency where the obtained approximations are not valid. For the 
band gap attributed to the shell resonance with index n = 1 the approximations based on 
the Foldy's equation predicts smaller gap width then that obtained with the self-consistent 
method. This is due to the fact that the lower limit is fixed on the shell resonance whereas 
in the self-consistent method it depends on the filling fraction. In Figure EJb) the shell 
radius is fixed and the lattice constant is varied. This results in the variation of the Brag 
frequencies. The obtained results are similar to those in Figure EJ^a). Again the derived 
approximations should accurately predict n = band gap width. 
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Figure 5: (Color online) Limits of the band gap n = and n = 1 obtained with the 
Foldy's approximations (|16p (' — '), MAE solution (|55|) ('• • •') and self-consistent method 
(fH7|) ('- • -'). The first Bragg frequency is plotted with solid line ('-') (a) Variation of the 
shell mid-surface radius a. (b) Variation of the lattice constant L. (c) Variation of the 
thickness 2h with filling fraction T ~ 0.4. (d) Variation of the Young's modulus E with 
filling fraction T ~ 0.4. 
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It can be observed in Figures [5j a) and (b), that n = band gap width quickly increases 
with the filling fraction. For T = 0.4 its width is approximately 300 Hz, while for T = 0.6 
it is 500 Hz. This trend is less pronounced for n = 1 band gap which remains quite narrow 
even in dense arrays. For example, the approximations based on the Foldy's equation 
estimate the band gap width as 10 Hz for T = 0.4 and 12 Hz for T = 0.6 whereas the 
self-consistent method gives 20 Hz and 35 Hz, respectively. 

In Figure [5jc) the dependence on the shell thickness is illustrated. The variation of 
the shell thickness also leads to the change in the filling fraction. However, due to the 
small increment in the outer radius of the elastic shell it is assumed that the filling fraction 
is fixed (for example, T ~ 0.4 for 2h = 2 mm and T ~ 0.37 for 2h = 0.25 mm). It 
is expected that the axisymmetric resonance and the corresponding band gap are shifted 
to the higher frequencies with decrease in the shell thickness. The approximations must 
become less accurate as the resonance approaches the Bragg frequency so that the results 
are analysed for 2h > 0.1 mm. The increase in the shell thickness results in the substantial 
low- frequency shift of n = band gap. It is predicted that the width of n = band gap 
observed around 1.1 kHz is approximately 250 Hz for the shell thickness 2h = 0.25 mm. For 
the shell thickness 2h = 1.5 mm the band gap is observed around 550 Hz with width 100 
Hz. It is also noted that the width of n = 1 band gap increases as the thickness decreases. 

The dependence of the upper and lower band gap limits on the elastic modulus of the 
shell is shown in Figure [5jd). The variation of Young's modulus (i.e. E £ [0.1, 20] MPa) 
results in a shift of the band gap limits. The approximations based on the Foldy's equation 
underestimate the upper limit of n = band gap. This is similar to the results shown in 
Figures [5ja), (b) and (c). It is observed that with decrease of the Young's modulus the 
band gaps appear at lower frequencies. In the given Young's modulus interval the width of 
n = band gap experiences little change (i.e. the width is approximately 250 Hz) whereas 
n = 1 band gap width becomes smaller as the Young's modulus decreases. For example, 
the approximations based on the Foldy's equation estimate the width of n = 1 band gap 
as 20 Hz for E = 10 MPa and 10 Hz for E = 1 MPa. The self-consistent method predicts 
the band gap width double of that obtained with the Foldy's equation. 

6 Conclusions 

In this paper the limits of the band gaps observed in a doubly periodic array of elastic 
shells and related to the scatterer's resonances are derived. These band gaps are observed 
below the first Bragg band gap. It is noted that the upper limit of the band gap attributed 
to the axisymmetric resonance n = is in the vicinity of the first Bragg band gap that 
restricts the use of the Foldy's equation. It is shown that this approach can only be 
used for the substantially small filling fractions (i.e. < 0.4) when the periodicity has 
little effect. The improved equation for the upper limit is derived by using the matched 
asymptotic expansions. This approach includes the periodic effects that results in the 
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accurate approximations for the densely packed periodic array (i.e. T > 0.4). In the 
low-frequency regime the alternative to the Foldy's equation is derived by using the self- 
consistent method. It is shown that these results can be used for the filling fractions 
bigger than 0.4. The methods described here can be easily adapted to derive the analytical 
expressions for low frequency band gap boundaries in arrays of circular resonant scatterers 
of different nature (for instance split rings or composite scatterers |22j). 
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